Trebern-Launay etal. BMC Medical Research Methodology 201 3, 13:102 
http://www.biomedcentral.eom/1 471-2288/1 3/1 02 



Medical Research Methodology 



RESEARCH ARTICLE Open Access 



Comparison of the risk factors effects between 
two populations: two alternative approaches 
illustrated by the analysis of first and second 
kidney transplant recipients 

Katy Trebern-Launay 1 ' 2 , Magali Giral 2 , Jacques Dantal 2 and Yohann Foucher 1 * 



Abstract 

Background: Whereas the prognosis of second kidney transplant recipients (STR) compared to the first ones has 
been frequently analyzed, no study has addressed the issue of comparing the risk factor effects on graft failure 
between both groups. 

Methods: Here, we propose two alternative strategies to study the heterogeneity of risk factors between two groups 
of patients: ft) a multiplicative-regression model for relative survival (MRS) and (ii) a stratified Cox model (SCM) 
specifying the graft rank as strata and assuming subvectors of the explicatives variables. These developments were 
motivated by the analysis of factors associated with time to graft failure (return-to-dialysis or patient death) in second 
kidney transplant recipients (STR) compared to the first ones. Estimation of the parameters was based on partial 
likelihood maximization. Monte-Carlo simulations associated with bootstrap re-sampling was performed to calculate 
the standard deviations for the MRS. 

Results: We demonstrate, for the first time in renal transplantation, that: (i) male donor gender is a specific risk factor 
for STR, (ii) the adverse effect of recipient age is enhanced for STR and (Hi) the graft failure risk related to donor age is 
attenuated for STR. 

Conclusion: While the traditional Cox model did not provide original results based on the renal transplantation 
literature, the proposed relative and stratified models revealed new findings that are useful for clinicians. These 
methodologies may be of interest in other medical fields when the principal objective is the comparison of risk factors 
between two populations. 



Background 

In patients facing a first allograft loss, repeat kidney 
transplantation provides a better chance for both long- 
term survival and quality of life than a return to dialysis 
[1,2]. The prognosis of second kidney transplant recipi- 
ents (STR) compared to first kidney transplant recipients 
(FTR) has been frequently studied. The older literature 
tends to conclude that STR have a worse prognosis than 
FTR [3,4]. However, recent analyses with adjustments 
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for confounding factors have challenged this generally 
accepted idea [5,6], with the exception of one study [7]. By 
modelling the time-dependent hazard between FTR and 
STR, we recently demonstrated that STR have a higher 
risk of graft failure than FTR, but this excess risk appears 
several years after transplantation [8]. According to this 
literature, one can accept that the excess risk for STR com- 
pared to FTR is negligible considering the improvements 
in life expectancy and quality of life compared to dialysis 
therapy. Nevertheless, as the demand for kidney trans- 
plants largely exceeds the supply [9], it is necessary to 
evaluate the differences in risk factors between STR and 
FTR so as to improve graft allocation. 
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For this purpose, traditional survival models can be used 
by merging STR and FTR. Nevertheless, one can notice 
two important limitations of these traditional approaches. 
Firstly, the comparisons of risk factors between both 
groups would imply testing all the interactions with the 
graft rank. Secondly, STR-specific explicative variables 
(survival time of the first transplant, transplantectomy 
or time in dialysis before re-transplantation) cannot be 
included despite the knowledge that their use would 
improve risk evaluation [5,6,10,11]. 

To overcome these limits, in this paper we propose 
two alternative strategies. The first one is an adaptation 
of a multiplicative-regression model for relative survival 
(MRS). This type of relative approach is often used to 
study the net survival of patients with cancer, i.e. the 
survival of patients if the only cause of death is related 
to the disease [12-17]. The principle of such additive- 
regression models is to introduce expected mortality rates 
by using life tables adjusted for gender, age and calendar 
year. Still using life tables, Andersen et al. [18] proposed 
a multiplicative-regression model. To our knowledge, the 
development of such methodology to endpoints other 
than mortality and with a reference group without a life 
time, has never been explored. 

Moreover, we propose a second method by adapting 
a stratified Cox model (SCM) specifying the graft rank 
as strata and assuming a vector of explicative variables 
decomposed into subvectors of variables that enter either 
in the reference hazard only, or in the relative hazard only, 
or in both groups but with common or separate effects. 

Methods 

Study population 

Second transplant recipients (STR) constituted the rel- 
ative group of interest. Recipients older than 18 years 
at the date of transplantation between 1996 and 2010 
were selected from the French DIVAT (Donnees Informa- 
tisees et VAlidees en Transplantation - www.divat.fr/en) 
multicentric prospective cohort [19]. Codes were used to 
assure donor and recipient anonymity and blind assay. 
The 'Comite National Informatique et Liberie' approved 
the study (N° CNIL 891735) and written informed con- 
sent was obtained from the participants. Only recipients 
with a maintenance therapy with calcineurin inhibitors, 
mammalian target of rapamycin inhibitors or belata- 
cept, in addition to mycophenolic acid and steroids, were 
included. Simultaneous transplantations were excluded. 
Recipients with at least one missing data for all the vari- 
ables taken into account in the expected hazard were 
excluded. The same criteria were applied to the reference 
group composed of first transplant recipients (FTR). The 
principal outcome was the time between transplantation 
and graft failure, which was the first event between return 
to dialysis and patient death with a functioning graft. 



The multiplicative-regression model for relative survival 
(MRS) 

Let the individuals be indexed by j (J = 1, n e ) in the ref- 
erence group and by i (i = 1, n r ) in the relative group. n e 
and n r represent the sample sizes of the respective groups. 
We note h°(fi\Zi) the observed instantaneous hazard func- 
tion at time t\ for the /th individual, where Zi is the vector 
of explicative variables. The observed hazard of the ith 
subject of the relative group can be decomposed in the 
multiplication of two hazards [18,20] 



h 0 (t i \z i ) = h\t i \z e i ) A'ft-K) 



(1) 



where h e (ti\zf>) is the expected hazard for an individual in 
the reference group with similar characteristics to the ith 
individual. z\ is a subset of Zi and represents these com- 
mon characteristics. h r (ti\z!) is the relative hazard with z\ 
being a subset of Z{. 

Parameters of the expected hazard can be estimated 
assuming a semi-parametric and proportional hazards 
(PH) model [21]. For the ;th individual (j = 1, n e ) 

h e (tj\zp=h e 0 (tj)exp(p e zp (2) 

where h e Q {tj) is an unknown expected baseline hazard 
function and fi e is the vector of regression coefficients 
associated with zfj. The estimations fi e are obtained 
by maximizing the partial log-likelihood among the n e 
patients of the reference group 



logP£ e G8 e ) = X> 



\ht k >ti / 



(3) 



where 8j equals 1 if the failure was observed for the yth 
subject and 0 otherwise. The estimation of the variance- 
covariance matrix, V(P e ), is obtained via the correspond- 
ing information matrix. 

Parameters of the relative hazard can also be estimated 
using a semi-parametric PH model. For the ith individual 
(i = 1, n r ), the instantaneous hazard is defined as 



h r (t i \z i ) = h r 0 (t i )exp(p r z t ;) 



(4) 



where h r 0 (ti) is an unknown relative baseline hazard func- 
tion and fi r is the vector of regression coefficients asso- 
ciated with z r t . By adapting the partial likelihood function 
(3) and assuming the previous estimations of expected 
parameters as constants, the regression coefficients fi r are 
estimated by maximizing 



\ogVC r {p r ) = J2 *i 



i=l 



-log J2 ex P(^4) 

\k:t k >ti ) 



(5) 
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Of note, for explicative variables not taken into account 
in the expected hazard in the model (1), exp(/3 r ) repre- 
sents the observed hazard ratios (HR) in the reference 
group, just as exp(^ r ) estimated from the PH model (4) 
among the relative group. In contrast, for explicative vari- 
ables taken into account in the expected hazard model in 
the model (1), exp(/3 r ) represents the weighting factors 
between the expected HR, i.e. exp(/3 e ), and the observed 
HR in the relative group, i.e. exp(j3 e ) x exp(/T). In other 
words, for explicative variables involved in both models 
(2) and (4), fi r = 0 means that the variable has the same 
effect in both groups. If /3 r > 0, the hazard ratio increases 
in the relative group compared to the reference group. If 
P r < 0, the hazard ratio decreases. 

Nevertheless, in contrast to the traditional relative sur- 
vival models based on life tables, the expected hazards 
cannot be reasonably assumed as constants since the cor- 
responding parameters were estimated from the reference 
sample. To take into account the variability associated 
with the expected model (2) in the estimation of the 
relative model (4), we used Monte-Carlo simulations asso- 
ciated with bootstrap resampling [22]. At each of the B 
iterations (b = 1, . . . ,B), this procedure can be divided 
into the following steps 

(a) Generation of a vector of parameters using the 
multivariate normal distribution Af (j$ e , V(/3 e )^ 
obtained from the maximisation of the partial 
likelihood (3). This first step takes into account the 
variance of the expected hazard. 

(b) Generation of a bootstrap sample from the relative 
sample comprising n r subjects. This second step 
takes into account the variance due to 
sample-to-sample fluctuation. 

(c) From this bootstrap sample, the model (4) is 
estimated by maximizing (5) in which the simulated 
parameters /?|* are used instead of fi e . fi r b is the 
resulting estimation of the relative regression 
coefficients. 

Means, standard deviations and 95% confidence inter- 
vals can be calculated from the B estimations of /3|*. 

The stratified Cox model (SCM) 

For the /th individual of the overall sample (i = 1, n), i.e. 
the reference and the relative group together (n - n e + n r ), 
let Zi be the vector of covariates that are applied to the 
model for either first or second transplant recipients. Let 
k be the indicator of the strata with k = e for the reference 
group and k = r for the relative group. The stratified Cox 
model is given by 



hkifi\Zi) = h k)0 (ti) exp(pZi); for k = e, r. (6) 



with hfcfiiti) the baseline hazard function in the strata /<. 
The vector zi can be decomposed into four different sub- 
vectors: (i) the vector z\ of explicative variables that enter 
in the reference hazard only (their values equal 0 if k = r); 
(ii) the vector z\ of explicative variables that enter in the 
relative hazard only (their values equal 0 if k — e); (Hi) 
the vector z\ of explicative variables associated to both 
groups; and (iv) the vector z? a subvector of explicative 
variables included in z\ but with separate effects. The 
model (6) can be developed as follows 

h(ti\Zi) = h k>0 (ti) expWi + e> c z\ + p s ^ ir + p r z\y, 
for k = e, r. 

(7) 

with /3 = (fi e , fi c , fi s , fi r ) the vector of regression coeffi- 
cients and 8i r equals 1 if the subject i belongs to the strata 
k = r and 0 otherwise. Then, for the reference group, we 
obtain 

h e (ti\Zi) = h e , 0 (ti) exp^ e z\ + P c z<l) (8) 
and for the relative group, we obtain 

h r (ti\Zi) = h r , 0 (ti) expWi + PA + P V Z\) (9) 

Therefore, exp(fi e ) represents the HR associated with 
specific variables for the reference group. And exp(/3 r ) 
represents the HR associated with specific variables for 
the relative group. For variables only included in z c (and 
not in z s ), exp(fi c ) represents the common HR in both 
groups associated with z c . For variables included in both 
z c and z s , the HR associated with z c equals exp(/3 c ) in the 
reference group and exp(/3 c + /3 s ) in the relative group. 
Thus, exp(/3 s ) represents the weighting factors between 
the expected HR, i.e. exp(fi c ), and the observed HR in the 
relative group, i.e. exp(fi c ) x exp(fi s ). 

Evaluation of the proportional hazards assumption 

In models (2), (4) and (7), hazard proportionality was 
checked for each explicative variable by plotting log- 
minus-log survival curves obtained by the Kaplan and 
Meier estimator [23] and by testing the scaled Schoenfeld 
residuals [24] separately in the reference and relative sam- 
ples. If the observed hazard ratios are constant regardless 
of time in both groups, the ratio between both observed 
hazard ratios, i.e. the weighting factor exp(/3 r ), will also be 
constant. 

Software 

All statistical analyses were performed using R version 
2.15.1 [25]. The proposed multiplicative-regression model 
for relative survival was implemented in an R package 
MRsurv available at www.divat.fr/en/softwares/mrsurv. 
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The adapted stratified Cox model was implemented by 
using the R package survival (function coxph, option 
strata). 

Results 

Description of the cohort 

641 STR potentially made up the relative group of inter- 
est, but 75 STR (11.7%) with missing data for explicative 
variables of the expected hazard were excluded. Finally, 
566 STR were included in the group of interest. The mean 
follow-up was 3.1 years with a maximum of 13.1 years. 
During the observation period, 72 returns to dialysis and 
34 deaths were observed. We identified 2462 FTR who 
met the inclusion criteria. We excluded 256 FTR (10.3%) 
with one missing data for at least one of the variables taken 
into account in the expected hazard model. Finally, 2206 
FTR made up the reference group. The mean follow-up 
was 3.4 years with a maximum of 13.7 years. During the 
observation period, 191 returns to dialysis and 109 deaths 
were observed. 

The demographic and baseline characteristics at the 
time of transplantation are presented in Table 1. Regard- 
less of the group, the majority of patients received a 
transplant from a deceased donor and the recipient gen- 
der was comparable between groups. However, STR were 
younger and their transplants were provided by younger 
donors. Recurrent nephropathy, past history of cardiac 
disease, hepatitis and malignancy were more frequent 
among STR, but STR had less diabetes and were less 
likely to be obese at the time of transplantation. Compared 
to FTR, STR received better HLA-matched transplants, 
but their cold ischemia time was longer and they were 
more immunized against HLA class I and class II antigens 
(historical Panel Reactive Antibodies) than FTR. They 
were also more frequently exposed to induction therapy 
with a lymphocyte-depleting agent. 

Among FTR meeting the inclusion criteria, some 
patients were also part of the STR group as they had 
received two transplants during the observation period. 
These 37 patients, who were included in both cohorts, 
represented 2% and 7% of the FTR and STR groups respec- 
tively. Given the large number of explicative variables, it 
seemed reasonable to assume conditional independence 
of the two transplantations of a given patient. In order 
to validate this assumption, we performed a frailty Cox 
model [26] based on the 37 individuals who were included 
in both groups. The frailty term was assumed to be 
Gamma distributed. The variance of the random vari- 
able was estimated at 5.10 -9 (p = 0.9948). Therefore, no 
intra-individual dependency was demonstrated. In order 
to validate the robustness of the results, we also performed 
both models after exclusion of the 37 STR also included 
in FTR. These results are presented in Additional files 1 
and 2. 



Analysis of risk factors in the FTR sample 

As previously illustrated in Table 1, it is well-established 
that FTR and STR are not intrinsically comparable. Thus, 
for the analysis of risk factors in the FTR population, 
adjustments were made (i) for all of the possible pre- 
or per-transplant immunological and non-immunological 
confounding factors according to experts and (ii) for all 
the baseline parameters differentially distributed between 
FTR and STR. All together, the expected hazard of graft 
failure was estimated according to recipient age and gen- 
der, causal nephropathy, comorbidities (including history 
of diabetes, hypertension, cardiac or vascular disease, 
dyslipemia, hepatitis B or C, and malignancy), obesi- 
ty, pre-transplant immunization (panel reactive antibody, 
PRA) against class I and class II antigens), donor age, 
deceased or living donor status, Epstein-Barr Virus (EBV) 
serology, period of transplantation, level of HLA-A-B-DR 
mismatches, induction therapy and cold ischemia time. 
This modelling is explained in detail in the paper by 
Trebern-Launay et al. [8] . The final multivariate model in 
the reference group of FTR is presented in the first three 
columns of Table 2. 

Relative hazard modelling in the STR group using the MRS 

A first selection of variables was performed (p < 0.20) fol- 
lowed by a step-by-step descending procedure (Wald test 
with p < 0.05). In line with the requirements of additive- 
regression models, adjustments were forced for recipient 
gender and age and transplantation period. All the vari- 
ables were categorized in order to avoid any log-linearity 
assumption and to obtain interpretable results. 

The final relative model is presented in the last 
three columns of Table 2. Expected HR previously esti- 
mated in FTR are presented in the first columns to 
enable a direct comparison between FTR (Cox model) 
and STR (relative model). Donor gender and wait- 
ing time before retransplantation were not taken into 
account in the expected hazards for FTR. Donor gen- 
der was not a significant risk factor for FTR and 
waiting time is by definition a specific factor for 
STR. More precisely, we estimated a 1.5-fold increase 
in risk of graft failure for STR with grafts from 
males compared to STR with grafts from females 
(p = 0.0320). Moreover, STR who waited more than 3 
years in dialysis before retransplantation had a 1.9-fold 
increased risk compared to STR with a shorter waiting 
time (p < 0.0001). 

In contrast, the effect of recipient age and donor age 
seemed significantly different between FTR and STR 
(p < 0.05). More precisely, if we assumed a similar effect 
of recipient age between both groups, the expected HR 
associated with recipient age > 55 years would be 1.39 
in the STR group, regarding the HR observed in the FTR 
group. In fact, the relative model showed that this HR 
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Table 1 Demographic characteristics at the date of transplantation for (i) whole cohort and (ii) FTR and STR separately; 
the last three rows of the table concern covariates specific for STR 

All (N = 2772) FTR (N = 2206) STR (N = 566) 

Demographic characteristics Number (%) Number (%) Number (%) p-value 



Transplantation period < 2005 


594 


(21.4) 


Recipient > 55 years of age 


1175 


(42.4) 


Male recipient 


1705 


(61.5) 


Recurrent causal nephropathy 


906 


(32.7) 


History of diabetes 


306 


(11.0) 


History of hypertension 


2263 


(81.6) 


History of vascular disease 


352 


(12.7) 


History of cardiac disease 


903 


(32.6) 


History of dyslipemia 


799 


(28.8) 


History of malignancy 


228 


(8.2) 


History of hepatitis B or C 


168 


(6.1) 


Recipient BMI > 30 kg.rrr 2 


263 


(9.5) 


Positive anti-class 1 PRA 


706 


(25.5) 


Positive anti-class II PRA 


733 


(26.4) 


Donor > 55 years of age 


1172 


(42.3) 


Deceased donor 


2470 


(89.1) 


Donor serum creatinine > 1 33 /xmol/l 


342 


(12.5) 


Positive donor EBV serology 


2613 


(94.3) 


HLA-A-B-DR incompatibilities > 4 


365 


(13.2) 


Cold ischemia time > 24h 


754 


(27.2) 


Lymphocyte-depleting induction 


1223 


(44.1) 



First graft survival < 1 year 

Waiting time before reg raft > 3 years 

Transplantectomy of the first graft 

p-values were obtained by using the Chi-square statistic. 

was 1.6-fold higher for STR compared to FTR (CI95% 
= [1.01-2.72], p = 0.0480). Similarly, the effect of donor 
age > 55 years was nearly two fold lower for STR than 
for FTR (CI95% = [0.33-0.99], p = 0.0440), while it was 
identified as a significant risk factor for FTR (HR = 
1.34, p = 0.0313). Of note, the relationship between the 
recipient gender and the risk of graft failure was not 
found to be significantly different between FTR and STR 
(p = 0.0720). 

Relative hazard modelling in the STR group using the SCM 

As an alternative, we performed the SCM based on the 
same variables as those used in the previous MRS. Donor 
gender and waiting time before retransplantation were 
included in variables applied only in the relative part, 
i.e. z r . The other four variables (transplantation period, 
recipient gender and age, and donor age) were included in 
z c and z s to evaluate the difference in their effect between 
both groups. The results are presented in Table 3. 



457 


(20.7) 


137 


(24.2) 


0.0806 


994 


(45.1) 


181 


(32.0) 


<0.0001 


1362 


(61.7) 


343 


(60.6) 


0.6536 


666 


(30.2) 


240 


(42.4) 


<0.0001 


269 


(12.2) 


37 


(6.5) 


0.0002 


1804 


(81.8) 


459 


(81.1) 


0.7545 


272 


(12.3) 


80 


(14.1) 


0.2804 


686 


(31.1) 


217 


(38.3) 


0.0012 


661 


(30.0) 


138 


(24.4) 


0.0104 


147 


(6.7) 


81 


(14.3) 


<0.0001 


96 


(4.4) 


72 


(12.7) 


<0.0001 


235 


(10.7) 


28 


(4.9) 


<0.0001 


355 


(16.1) 


351 


(62.0) 


<0.0001 


319 


(14.5) 


414 


(73.1) 


<0.0001 


973 


(44.1) 


199 


(35.2) 


0.0002 


1940 


(87.9) 


530 


(93.6) 


0.0002 


279 


(12.8) 


63 


(11.4) 


0.3807 


2087 


(94.6) 


526 


(92.9) 


0.1540 


326 


(14.8) 


39 


(6.9) 


<0.0001 


552 


(25.0) 


202 


(35.7) 


<0.0001 


793 


(35.9) 


430 


(76.0) 


<0.0001 






131 


(24.1) 








272 


(49.8) 








220 


(38.9) 





Estimations and corresponding 95% confidence inter- 
vals were very similar to those obtained in the MRS. 
Indeed, as in the previous model, we estimated a 2-fold 
increase in risk of graft failure for STR who waited more 
than 3 years in dialysis before retransplantation compared 
to STR who waited less than 3 years (p = 0.0019). The 
relationship between the donor gender and the risk of 
graft failure among STR was similar to that obtained in 
the MRS but was not found to be significant (HR = 1.51, 
p = 0.0674). 

Transplantation period, recipient gender, recipient age 
and donor age were included in variables applied in both 
models. Results were also concordant with the MRS. 
For the four explicative variables, estimations and corre- 
sponding 95% confidence intervals were similar to those 
obtained in the MRS. However, conversely to the MRS, 
recipient age and donor age were not found to be signif- 
icantly differently associated with the risk of graft failure 
between the two groups. 
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Table 2 Multivariate Cox model for FTR and results of the MRS in the STR group : (i) the first three columns provides the 
results of the multivariate Cox model analysis of graft failure risk factors for FTR (N = 2206); (ii) the next three columns 
provide the results of the relative survival model based on 540 STR (26 recipients presenting missing data for the waiting 
time before re-transplantation were excluded) 






Cox model in 
the FTR group 




MRS in 
the STR group 




HR 


95% CI 


p-value 


HR 95% CI p-value 


Variables entering in the model for FTR only 










Causal nephropathy (recurrent/ non recurrent) 


1 .24 


0.96-1 .59 


0.0987 




History of diabetes (positive/ negative) 


1 .34 


0.96-1 .85 


0.0819 




HKtnrv nf hvnprtpn^ion fnn^itivp / npnativpl 


0.77 


0.57-1 .05 


0.0986 




History of cardiac disease (positive / negative) 


1.41 


1.1 1-1.79 


0.0051 




HKtnrv nf \/rKn ilsr Hkps^p fnrKitivp / npnsth/pl 


1 .1 0 


0.81-1 .51 


0.5351 




History of dyslipemia (positive / negative) 


1.12 


0.87-1.45 


0.3828 




History of hepatitis B/C (positive / negative) 


0.82 


0.45-1 .47 


0.4969 




History of malignancy (positive / negative) 


1.25 


0.84-1.86 


0.2698 




Body mass index (> 30 kg.m-2 / < 30 kg.m-2) 


1.58 


1.12-2.14 


0.0084 




Anti-class 1 PRA (positive / negative) 


1.45 


1.07-1.97 


0.0182 




Anti-class II PRA (positive / negative) 


1.09 


0.78-1.52 


0.6299 




Donor status (deceased/living) 


2.50 


1.41-4.43 


0.0016 




Donor EBV serology (positive / negative) 


1.65 


0.98-2.78 


0.0606 




Number of HLA-A-B-DR mismatches (> 4 / < 4) 


1.30 


0.97-1.76 


0.0824 




Induction therapy (depleting / non depleting) 


0.79 


0.60-1.05 


0.1091 




Cold ischemia time (> 24 h / < 24 h) 


1.29 


1.01-1.66 


0.0441 




Variables entering in both models 










Transplantation period (< 2005 / > 2005) 


1.33 


0.97-1.82 


0.0693 


0.97 0.55-1.74 0.9360 


Recipient gender (male / female) 


1.17 


0.91-1.51 


0.2186 


0.61 0.38-1.05 0.0720 


Recipient age (> 55 years / < 55 years) 


1.39 


1.05-1.83 


0.0204 


1.65 1.01-2.72 0.0480 


Donor age (>55 years / <55 years) 


1.34 


1.03-1.74 


0.0313 


0.59 0.33-0.99 0.0440 



Variables entering in the model for STR only 



Donor gender (male /female) - - - 1.53 1.03-2.48 0.0320 

Waiting time before reg raft > 3 years - - - 1.92 1.22-3.00 <0.0001 

PRA, panel reactive antibody; EBV, Epstein-Barr virus; HLA, human leukocyte antigen. 



Discussion 

Although the comparison of survival between first and 
second kidney transplants has been frequently performed, 
no study has addressed the issue of comparing the risk fac- 
tors associated with the time to graft failure between both 
groups. Understanding the factors influencing the long- 
term evolution of STR compared to FTR would benefit the 
medical management of graft attribution by identifying 
patients with the best chances. 

The absence of literature focusing on this question 
may be partially explained by the methodological issues 
associated with such studies. Indeed, the Cox model is 
classically used to explore risk factors influencing graft 
survival and interactions can be included to evaluate risk 
factor differences between FTR and STR. However, this 



approach has several limitations. Firstly, it implies testing 
interactions between the graft rank and each explicative 
variable, increasing the number of parameters and mak- 
ing interpretations difficult. Secondly and certainly more 
importantly, only covariates common to both groups can 
be taken into account. This excludes explicative variables 
specific for one group. Concerning our application, this 
constitutes a limitation as several STR-specific explicative 
variables are known to be associated with second graft 
prognosis: the first graft transplantectomy [10], the first 
graft survival duration [6,11] or the time in dialysis before 
re-transplantation [5]. 

This paper describes two alternative models to 
overcome these difficulties. Firstly, the adaption of a 
multiplicative-regression model for relative survival 
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Table 3 Results of the stratified Cox model based on 2746 patients with 2206 FTR and 540 STR (26 STR presenting 
missing data for the waiting time before re-transplantation were excluded) 

FTR strata STR strata 

HR 95% CI p-value HR 95% CI p-value 

Variables entering in z e only 

Causal nephropathy (recurrent / non recurrent) 

History of diabetes (positive / negative) 

History of hypertension (positive / negative) 

History of cardiac disease (positive / negative) 

History of vascular disease (positive / negative) 

History of dyslipemia (positive / negative) 

History of hepatitis B/C (positive / negative) 

History of malignancy (positive / negative) 

Body mass index (> 30 kg.m-2 / < 30 kg.m-2) 

Anti-class I PRA (positive / negative) 

Anti-class II PRA (positive / negative) 

Donor status (deceased/living) 

Donor EBV serology (positive / negative) 

Number of HLA-A-B-DR mismatches (> 4 / < 4) 

Induction therapy (depleting / non depleting) 

Cold ischemia time (> 24 h / < 24 h) 
Variables entering in z c and z s 

Transplantation period (< 2005 / > 2005) 

Recipient gender (male / female) 

Recipient age (> 55 years / < 55 years) 

Donor age (>55 years / <55 years) 
Variables entering in z r only 

Donor gender (male /female) - - 1.51 0.97-2.36 0.0674 

Waiting time before regraft > 3 years - - 1.99 1.29-3.07 0.0019 

PRA, panel reactive antibody; EBV, Epstein-Barr virus; HLA, human leukocyte antigen. 



1.12 


0.90-1.39 


0.3031 








1.28 


0.95-1.72 


0.1001 


- 


- 


- 


0.88 


0.68-1.14 


0.3270 


- 


- 


- 


1.35 


1.10-1.66 


0.0042 


- 


- 


- 


1.11 


0.85-1.47 


0.4433 


_ 


_ 


_ 


1.19 


0.95-1.49 


0.1263 








0.97 


0.65-1.46 


0.9091 








1.21 


0.87-1.67 


0.2646 








1.54 


1.13-2.08 


0.0057 








1.39 


1.07-1.82 


0.0153 








0.98 


0.73-1.31 


0.8857 








2.06 


1.27-3.36 


0.0036 








1.65 


1.07-2.54 


0.0235 








1.33 


1.01-1.75 


0.0397 








0.88 


0.70-1.11 


0.2742 








1.20 


0.97-1.49 


0.0894 









1.42 1.09-1.86 0.0099 

1.17 0.91-1.50 0.2200 

1.36 1.03-1.78 0.0274 

1.36 1.04-1.77 0.0238 



0.94 0.54-1.64 0.8295 

0.63 0.40-1.02 0.0581 

1.60 0.95-2.72 0.0785 

0.60 0.35-1.05 0.0725 



allows a direct comparison of risk factors between two 
groups of patients without presupposing the role of each 
variable, i.e. common, different or specific relationships. 
The corresponding semi-parametric models are the Cox 
model for the expected hazard and the multiplicative - 
regression for the relative hazard. The main difficulty 
and limit of these models is the estimation of standard 
deviations which were obtained by Monte-Carlo sim- 
ulations associated with bootstrap re-sampling. In this 
multiplicative modelling, the regression coefficients are 
straightforward to interpret in terms of their interactions. 
We propose an R package for a simple way of using the 
model. 

Secondly, we demonstrated that a stratified Cox model 
specifying the graft rank as strata may be fitted to 
take into account STR-specific variables as a subvector 
of variables that enter in the model for STR only. In 



addition, some variables would enter either in the model 
for FTR only or in both models (with common or sep- 
arate effects). The main limit to this approach is that 
the corresponding structure presupposes knowledge of 
variables potentially applicable to both models (in con- 
trast with the relative model) unless testing a very large 
number of models. Indeed, whereas explicative variables 
entering in a single model (for FTR or STR) would 
easily be clinically assumed, those applicable to both 
models and with common or separate effects are not 
known in advance. Nevertheless, the SCM can be sim- 
ply estimated by maximising a single partial likelihood 
function. 

As expected, the results were concordant between both 
approaches. Regression coefficients were similar while 
standard deviations appeared a little smaller with the MRS 
approach. The results showed that male donor gender 
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and a long waiting time before retransplantation were two 
specifioSTR risk factors: donor gender was not signifi- 
cantly associated with the risk of graft failure in the FTR 
population and the waiting time before retransplantation 
was only related to STR by definition. The interpretations 
were similar to hazard ratios from a Cox model preformed 
on the STR group. 

Conversely, two explicative variables appeared to be dif- 
ferently associated with the risk of graft failure between 
STR and FTR. More precisely, we showed for the first 
time that the adverse effect of recipient age was enhanced 
for STR as compared to FTR. The main clinical expla- 
nation is a cumulative effect of the risk factors for 
STR, in particular because of the cumulative exposure 
to immunosuppressive drugs during the first transplan- 
tation period. From a clinical point of view, this result 
may imply that clinicians should pay particular atten- 
tion to recipient age in second kidney transplantations. 
Also, for the first time to our knowledge, this study iden- 
tified an attenuation of the risk factor related to older 
transplants for STR as compared to FTR. Two explana- 
tions are: (i) an indication bias with only high-quality 
donors (without diabetes, hypertension or cardiovascular 
disease) proposed to STR; (ii) a higher non-HLA immu- 
nization in STR, explaining why graft failure is generally 
due to immunological phenomena rather than transplant 
quality. 

Although we illustrated the advantages of both alter- 
native approaches in renal transplantation, this method- 
ology may be useful in number of other clinical and 
epidemiological applications. For practical use, we pro- 
pose an R package to compute the MRS. The adaptation 
of the SCM can be computed by using many statisti- 
cal software. Of course, the aim of such models is not 
to replace traditional survival models, but rather to pro- 
vide a more suitable alternative when the main objective 
is to compare risk factors between two populations, in 
particular when population-specific covariates need to be 
included. 

As always, there are several avenues worth exploring 
from this work. First, both models can be generalized for 
time-dependent explicative variables by adapting the like- 
lihood functions as proposed by Therneau and Grambsch 
[27] (chapter 5, pages 111-115). Second, both models 
assumed the independence of FTR and STR. While this 
assumption was evaluated by using a frailty model among 
the 37 individuals common in both groups, a low statisti- 
cal power may explain the non-rejection of this indepen- 
dence hypothesis. To ensure the validity of our results, 
we reperformed both models after exclusion of the cor- 
responding 37 STR also included in FTR. The results 
presented in additional files showed the robustness of 
the results. Third, other strategies for variable selection 
can be adapted, such as partial likelihood generalization. 



Finally, further work is needed to develop a Goodness-of- 
fit statistic for the MRS approach, in particular concerning 
the proportional hazards assumption. 

Conclusions 

MRS and SCM consitute two original approaches to com- 
pare risk factors between two populations. The advantage 
of MRS is to allow a direct modelling strategy but it is 
not straightforward to estimate the standard deviations. 
In contrast, SCM allows an overall estimation of param- 
eters and standard deviations but its structure presup- 
poses knowledge of the role of each explicative variable. 
This study also highlighted novel risk factor differences 
between first and second kidney transplant recipients. 
These results could help improve the management of 
patients waiting for a second graft. They may also encour- 
age the widespread use of this original methodology in 
other medical fields. 

Additional files 



Additional file 1 : Multivariate Cox model for FTR (N = 2206) and 
results of the MRS in the STR group after exclusion of the 37 STR also 
included in FTR, based on 507 STR (22 recipients presenting missing 
data for the waiting time before re-transplantation were excluded). 

Additional file 2: Results of the stratified Cox model after exclusion 
of the 37 STR also included in FTR, based on 271 3 patients with 2206 
FTR and 507 STR (22 STR presenting missing data for the waiting time 
before re-transplantation were excluded). 
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